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Abstract 

In this paper we derive some new and practical results on testing and interval estimation 
problems for the population eigenvalues of a Wishart matrix based on the asymptotic theory 
for block-wise infinite dispersion of the population eigenvalues. This new type of asymptotic 
theory has been developed by the present authors in Takemura and Sheena (2005) and Sheena 
and Takemura (2007a, b) and in these papers it was applied to point estimation problem of 
population covariance matrix in a decision theoretic framework. In this paper we apply it to 
some testing and interval estimation problems. We show that the approximation based on this 
type of asymptotics is generally much better than the traditional large-sample asymptotics 
for the problems. 

Key words and phrases: eigenvalues of covariance matrix, Wishart distribution, test on eigenvalues, 
interval estimation of eigenvalues 

1 Introduction 

Let S = (sij) be distributed according to Wishart distribution W p (n, £), where p is the dimension, 
n is the degrees of freedom and S is the covariance matrix. Let Ai > . . . > X p > denote 
the eigenvalues of X. In this paper we consider some testing and interval estimation problems 
for the eigenvalues of S. Our aim is to give practical solutions to the problems based on the 
asymptotic theory for block-wise infinite dispersion of the population eigenvalues. In view of 
the intractability of the finite sample exact distribution of sample eigenvalues, usually the large 
sample asymptotic approximation is used. There exists an extensive literature on improving the 
first-order large sample approximation by an asymptotic expansion (see Siotani et al. (1985) for 
a comprehensive treatment). However for a moderate or small value of the sample size n, the 
large sample asymptotic theory often gives a poor approximation. In these cases asymptotic 
expansions tend to give an even larger error. On the other hand, we find that approximation 
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based on the infinite dispersion of the population eigenvalues is more robust and gives a reasonable 
approximation even for a small sample size n. 

The first problem we consider in this paper is testing the one-sided null hypothesis on the mth 
population eigenvalue 



For testing it is natural to consider a one-sided rejection region based on the mth sample 

eigenvalue of l m of S. We show that the least favorable distribution is given by = Ai = • • • = A m 
and = A m+ i = • • • = \. This is exactly the situation covered by the asymptotic theory for 
block-wise infinite dispersion. Therefore it gives an explicit solution to the testing problem of 

rr{m) 
n ■ 

The second problem is the interval estimation for the largest population eigenvalue Ai in terms 
of the largest sample eigenvalue li of S. We will show that confidence interval based block- 
wise infinite dispersion gives much better coverage probability than the conventional large sample 
asymptotics. 

The third problem is testing the hypothesis of equality of the several smallest eigenvalues: 
A m +i = • • • = A p . This problem is important in determining the rank of the systematic part in a 
multivariate variance component model. We consider approximation to the null distribution of the 
likelihood ratio criterion under the block-wise infinite dispersion of population eigenvalues. Again 
this type of asymptotics gives much better approximation than the large sample asymptotics. 

The organization of the paper is as follows. In Section [2] we set up notations for the paper 
and give some preliminary results on the asymptotic theory for block-wise infinite dispersion of 
the population eigenvalues. In Section [3] we study the above three problems, 1) One-sided test for 
a population eigenvalue in Section 3.1; 2) Interval estimation for extreme eigenvalues in Section 
3.2; 3) Testing equality of the smallest eigenvalues in Section 3.3. 

2 Asymptotic Distribution of Normalized Sample Eigen- 
values 

Denote the spectral decompositions of I] and S by 



where G, T e 0(p), the group of p x p orthogonal matrices, and A = diag(Ai, . . . , A p ), L = 
diag(/i, . . . ,l p ) are diagonal matrices with the eigenvalues Ai > . . . > X p > 0, l\ > . . . > l p > of 
S and S, respectively. We use the notations A = (Ai, . . . , A p ) and I = (7i, . . . , l p ) hereafter. We 
guarantee the uniqueness (almost surely) of the decomposition (jSJ) by requiring that 



has positive diagonal elements. 

In Takemura and Sheena (2005) we considered what happens to appropriately normalized 
components of S if the population eigenvalues become infinitely dispersed, i.e., 




(1) 



S = TAT' 
S = GLG', 



(2) 
(3) 



G = (g tl ) = T'G 



(4) 



(A 2 /A 1 ,A 3 /A 2 ,...,A P /A p _i) ^0. 
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In Sheena and Takemura (2007b) we generalized the asymptotic result of Takemura and Sheena 
(2005) to the case when the population eigenvalues are block-wise infinitely dispersed. 
Let the population eigenvalues be parameterized as follows; 

j &a, if i = l,..., m, , v 

Ai ifz = m + l,...,p, l ° j 

where £j's are fixed and "asymptotic parameter" a and j3 vary. When we say population eigenvalues 
are "(two-)block-wise infinitely dispersed" , it means that 

P/a -> 0. (6) 

The above notation is used as a general notation including specific convergences (divergences) 
such as (a, (3) — > (oo, 1), (a, P) — > (1,0) and so on. More precisely, the operation lim^/ Q ^ /(a, 0) 
means lim^oo f(a iy Pi) with any specific sequences oti, Pi, % = 1,2,..., such that Pi /on — > as 
z — ► oo. 

As Sheena and Takemura (2007b) indicates, appropriate normalization for the sample eigen- 
values is given by 

U/a, in = l,.-.,m, 

\ Zi//3, if z = m + l,...,p, v ; 

while 01]) also serves as an appropriate normalization for sample eigenvectors. 

For the normalized population or sample eigenvalues, we use the following notations; 





— {til - ■ ■■, im)i 


6 


(£m+l , • • • , £p) 


, 




= (dx, . . . , d m ), 


d 2 


(^m+1, • • • , ^p, 


1, 


Si 


= diag(G, . . 




= diag(£ m+ i, . . . 


, 6p): 


£>i 


= diag(di, ...,d m ), 


£>2 


= diag(d m+ i, . . . 


,,d p ) 



Now we state the basic theorem on the asymptotic distributions of dx,d 2 . 
Theorem 1 Suppose that we have two independent Wishart distributions 

Wn ~ W m (n, Hi), W 22 ~ W p - m (n - m, S 2 ) 
and that their spectral decompositions are given by 

W lx = G n DxG' n , Dx = diag(di, . . . , d m ), dx = (dx,..., d mh 

W 22 = G 22 D 2 G' 22 , D 2 = diag(cT m +i, • • • , dp), d 2 = (d m+1 , ■■■,d p ), 

where Gxx £ 0(m), G 22 G 0(p — m), dx > • • • > d m , d m+1 > • • • > d p . Then as P/a — > 0, 

d~, z = 1,2. 

Proof. Using Lemma 1 of Sheena and Takemura (2007b), we prove the convergence of the 
moment generating function. Let 

(m p \ / p 

a' 
i=l i=m+l / \i=l 



where \9i\ < 3 1 min J ^ 1 , Vi Notice that (19) in Lemma 1 of Sheena and Takemura (2007b) is 
satisfied since 



x 



(TG,l,\,a,P)<expla- 1 J2h\0 i \+p- 1 ^ k\9i 

\ i=l i=m+l 
( m 1 

< exp Z~ l a- 1 hii 1 + Z- 1 ?- 1 iiST 1 



i=l 



i=m+l 



v, i=l / 

< exp (trS^GXG'A' 1 ) , VG G 0(p), VZ G {Z|/i > • • • > l p > 0}. 



For the last inequality, see e.g. Marshall and Olkin (1979) Ch.20.A.l. Since x(d, q, £, a, (3; T, H^) 
exp (X)i=i trivially we have 

x T (H^G(q U: q 22: 0), cZ, Q 21 , £) = exp (j^i) ■ 



,i=i 



Therefore we have 



lim 

/3/a-»0 



i=l 



E 



,i=i 



i=m+l 



3 Inference on Population Eigenvalues 

The asymptotic result in the previous section has possibly various applications for inference on 
the population eigenvalues. We give three inference problems as interesting applications. 

3.1 One-sided Test for Population Eigenvalue 

Consider the null hypothesis on the mth (m — 1, . . . ,p) population eigenvalue 



-«0 • A m d. A m 



against the alternative : \ m < \* m . Need for testing H arises in some practical cases, for 

example: 



In principal component analysis, A*(= A* = • • • = A*) may be a cut-off value and a test for 

Hq™* is repeatedly carried out starting from m — 1 until is rejected. This is one of the 
methods for deciding the dimension of the principal components. 
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• Let Xi(i = 1, . . . ,p) be the return of the ith asset in finance and x = (xi, . . . , x p ) is dis- 
tributed as the p-dimensional normal distribution N p (0, E). is equivalent to the asser- 
tion a'Ha < A*, Va = (a 1; . . . , a p ) such that ||a|| = 1. If is rejected, then it means that 
the group of assets x is stable in view of volatility since any portfolio among the group is 
never beyond A^ in its variance. 

A natural rejection region in testing is given by l m < ^(7) for a given significance level 7. 
The following lemma and Theorem [1] give the critical point ^(7). 

Lemma 1 For any positive c 

sup P A (l m < c) = lim Px{l m < c), 

(m) P^O 

where A = diag(Ai, . . . , X p ), Ai = ■ • • = A TO = A^, X m +i = • • • = A p = (3. 

Proof. According to Theorem 1 of Anderson and Das Gupta (1964), P\(l m < c) is a monoton- 
ically decreasing function with respect to each Aj, (i = 1, . . . ,p), hence 

P\(Jm < C) < P- A {lm < C), 

where [3 = X p . Furthermore P\{l m < c) is monotonically increasing as (3 goes to zero. 1 

Because of the result of Theorem [1] with a — 1, & = A^, (i — 1, . . . , m) and = 1, (i — 
m+l,...,p), 

limP A (Z m < c) =P(Z m < c), 

where l m is distributed as the smallest eigenvalues of W m (n, A^J m ). Therefore we have the 
following result. 

Theorem 2 Por testing hypothesis H^ 1 against H[ m \ a test with significance level 7 is given 
with the rejection region 

i m < i* m (r/), 

where ^(7) is the lower 100y% point of the smallest eigenvalue ofW m (n, X* m I m ). 

For analytic calculation of ^(7), see Thompson (1962), Hanumara and Thompson (1968). In the 
case m — 1, which is practically the most important, it is given by A*Xn(7); where Xnil) is the 
lower 1007% point of the x 2 distribution with the degree of freedom n. 

3.2 Interval Estimation of Extreme Eigenvalues 

In this subsection we present a new way of constructing a confidence interval for the extreme 
population eigenvalues. Let Ai < fi(l) be a one-sided estimated interval with confidence level 7. 
For example, in the second case in Section I3.1[ the maximum volatility in all possible portfolio 
among the assets x is estimated to be less than or equal to 

However if we use the exact finite distribution theory, it is not easy to find an appropriate 
fi(l) under a given 7 even if we only consider an interval of the simplest form Ai < C\li with some 
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constant C\. (Note that k/Xi, (i = 1, . . . ,p) is bounded in probability. See Lemma 1 of Takemura 
and Sheena (2005).) Therefore usually a large sample approximation is employed (e.g. Theorem 
13.5.1. of Anderson (2003)): 

v^-A.) -iV(0,2A, 2 ), i = l,...,p. 
Let z 1 denote the upper IOO7 percentile of the standard normal distribution. Since 



= P(li> (V2nz y + n)X 

we have an approximate confidence interval 

A x < (v2nz 7 + n)^ 1 ^, 



7 as n — > 00, 



with confidence level close to 7 for sufficiently large n. 

Now we propose an alternative approximation. Suppose m — 1 in Theorem [H then as (3 /a 
goes to zero, 

d\ — > di = W11. 

Since W n /£i ~ X 2 { n ), 

h d\ d 2/ \ 
T- = 7- X (n). 
^1 si 

as /3/a goes to zero. From this asymptotics, we can make an approximate interval 

Ai < (Xy(n)r%, (9) 

where x^{ n ) is the upper IOO7 percentile of \ 2 distribution with the degree of freedom n. The 
interval (Q has approximately 7 confidence level when /3/a is sufficiently close to zero. 

We are interested in how large n for (jSj) or how small for (jUJ) is required to get practi- 
cally sufficient approximations. Because of difficulty in theoretical evaluations, we carried out a 
simulation study with the fixed parameters p = 3, m — 1, £1 = £2 = ^3 = 1, a — 1, while we 
select different n's (5, 10, 20, 50, 100, 500, 1000) and (3's (1.0, 0.9, 0.8, 0.6, 0.5, 0.3, 0.1, 0.01, 
0.001). For each case, 50000 Wishart random matrices are generated. We present the results in 
Table 1. The numbers under L1(U1) indicate the ratio of the largest sample eigenvalue which fall 
within the interval (jSj) with 7 = 0.95(0.05), while those under U2(L2) show the similar ratio with 
respect to the interval ([9]). Numbers in bold indicate that they are within ±0.01 deviation from 
the desired value, hence the approximation may be good enough for many practical purposes. We 
can summarize the result as follows; 

1. In every case gives better approximation than (JHJ). 

2. Since (3 is as large as 0.3, (El) already gives a good approximation. In that sense, the 
approximated interval ([HD seems robust. When (3 is smaller or equal to 0.1, (El) works well 
even with small samples such as n = 5 or 10, while needs samples as large as 100 or 500. 
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Table 1: Approximated Interval Estimation 



p 


1 


0.9 


0.8 




Ul 


U2 


LI 


L2 


Ul 


U2 


LI 


L2 


Ul 


U2 


LI 


L2 


n = 5 


.402 


.322 


1.00 


1.00 


.341 


.267 


1.00 


1.00 


.276 


.209 


1.00 


1.00 


n= 10 


.416 


.345 


1.00 


1.00 


.327 


.264 


1.00 


1.00 


.247 


.192 


1.00 


1.00 


n = 20 


.419 


.361 


1.00 


1.00 


.301 


.250 


1.00 


1.00 


.207 


.165 


1.00 


1.00 


n = 50 


.416 


.373 


1.00 


1.00 


.256 


.222 


1.00 


1.00 


.153 


.129 


1.00 


1.00 


n = 100 


.413 


.380 


1.00 


1.00 


.212 


.189 


1.00 


1.00 


.124 


.109 


.999 


.999 


n = 500 


.405 


.389 


1.00 


1.00 


.118 


.111 


.999 


.999 


.080 


.074 


.984 


.981 


n = 1000 


.407 


.395 


1.00 


1.00 


.097 


.093 


.995 


.994 


.070 


.067 


.972 


.970 



P 


.6 


.5 


.3 




Ul 


U2 


LI 


L2 


Ul 


U2 


LI 


L2 


Ul 


U2 


LI 


L2 


n = 5 


.167 


.120 


1.00 


1.00 


.133 


.094 


1.00 


1.00 


.095 


.067 


1.00 


.998 


n= 10 


.139 


.104 


1.00 


1.00 


.111 


.083 


1.00 


.999 


.084 


.063 


.999 


.989 


n = 20 


.112 


.089 


1.00 


.998 


.094 


.073 


.999 


.995 


.075 


.058 


.990 


.974 


n = 50 


.089 


.075 


.995 


.991 


.080 


.068 


.988 


.980 


.069 


.058 


.974 


.961 


n = 100 


.078 


.068 


.984 


.978 


.072 


.064 


.975 


.968 


.062 


.055 


.965 


.957 


n = 500 


.064 


.059 


.963 


.959 


.059 


.055 


.959 


.955 


.057 


.053 


.958 


.954 


n = 1000 


.060 


.057 


.960 


.958 


.058 


.055 


.956 


.953 


.055 


.053 


.956 


.954 




P 


.1 


.01 


.001 




Ul 


U2 


LI 


L2 


Ul 


U2 


LI 


L2 


Ul 


U2 


LI 


L2 


n = 5 


.075 


.055 


1.00 


.975 


.068 


.049 


1.00 


.953 


.071 


.052 


1.00 


.951 


n = 10 


.072 


.055 


.992 


.959 


.068 


.052 


.989 


.951 


.065 


.048 


.988 


.950 


n = 20 


.066 


.053 


.979 


.956 


.066 


.052 


.976 


.951 


.063 


.050 


.974 


.948 


n = 50 


.061 


.051 


.966 


.953 


.058 


.048 


.964 


.950 


.059 


.050 


.965 


.950 


n = 100 


.057 


.050 


.961 


.951 


.057 


.051 


.961 


.951 


.056 


.049 


.960 


.951 


n = 500 


.055 


.051 


.954 


.950 


.054 


.050 


.956 


.952 


.053 


.049 


.954 


.950 


n = 1000 


.052 


.050 


.954 


.951 


.053 


.050 


.953 


.949 


.052 


.049 


.954 


.951 



3. When (3 is from 0.5 to 0.6, both approximations need large samples such as 500 or 1000. If 
(3 is larger than 0.6, they need samples larger than 1000 for a good approximation. 

Similarly we can make an interval estimation for the smallest eigenvalue, X p . Let m = p — 1 in 
Theorem [TJ, then ft/a goes to zero, 

dp — > d p = W 22 - 

Since W-nHv ~ X 2 ( n ~ P + 1), 

lp dp d 2 1 , i \ 

f- = — -> x (n-p + 1). 

C,p 

as P/ct goes to zero. Using this fact, we can estimate X p to lie in the interval 

A P < (4(n-p + l))-%, (10) 

7 



at approximately 7 confidence level when /3/cx is sufficiently close to zero. 

We now compare (jHJ) and (Q more closely in view of the known results on asymptotic expansion 
of distribution of sample eigenvalues. Let 

/w fljjX 

2 V^T 

The asymptotic expansion of A n up to the order n 1 / 2 is given by (see Sugiura (1973)) 

+ o(n~ 1 / 2 ). (12) 



^-1). (id 



Now suppose x i? z = 1, . . . , n, are independently and identically distributed as x 2 (l) distribution. 
The normalized variable 

1 

7=(Xj 1), i I, ... , 71 

has zero mean and unit variance. Let 

n I — \ 
. — , / n I \ . t ■ \ 

(13) 



The asymptotic expansion of B n up to the order n 1//2 is given by 

F Bn (t) = m-^(t 2 -l) + o(n-V 2 ). (14) 

Comparing (1121) and (1141) . we notice that if t 2 > 1, then the absolute value of the second term 
in (j!4p is smaller than that of f|T2|) by the margin 



i=2 



Since /1/A1 is asymptotically distributed as x 2 ( n ) when the largest population eigenvalue is in- 
finitely deviated from the others, h/Xi in (jTTj) is similarly distributed as J^Xj in (fTBl) . In this case 
(TH)]) vanishes and both expansions (1T21 . f[T4|) coincide. It is naturally conjectured that when the 
largest population eigenvalue Ai is positioned far away from the smaller eigenvalues, we can make 
an "easier" inference on Ai. The fact that the term f|T5|) shrinks in that situation supports this 
conjecture as well as our simulation results. 

3.3 Testing Equality of the Smallest Eigenvalues 

As in the introduction of Section 11.7.3 of Anderson (2003), the equality of the p — m smallest 
population eigenvalues 

A m +i = • • ■ = A p (say a 2 ) (16) 
is equivalent to the covariance structure 

S = $ + <T 2 I p , 



s 



where a positive semidefinite matrix with rank m, represents the variance-covariance matrix 
of a systematic part and u 2 I p arises from measurement error. If hypothesis (fTSl) is accepted, then 
it suggests that the systematic part might consist of m independent factors. Need for testing (ITS"]) 
also arises in principal component analysis when the dimension of principal components has to be 
decided. Once it is accepted and a 2 is sufficiently small, which might require another hypothesis 
testing, we could ignore the last p — m principal components. 

The likelihood ratio statistic for testing ({TBI) is given (see e.g. Theorem 9.6.1 of Muirhead 
(1982)) by 



=m+l l » 



and the critical region is V < 0(7) for a given significance level 7 . 

In order to give the critical point 0(7) , we traditionally make use of the asymptotic convergence 

— nlog V -—>■ x 2 {{p ~ m + 2)(p — m — l)/2), as n —>■ 00. (17) 

Bartlett adjustment and further refinement on the asymptotic result are found in Section 9.6 of 
Muirhead (1982). From this convergence, the approximate critical point is given as 



0(7) = exp^— n 1 Xy((p — m + 2)(p — m — l)/2] 



On the other hand we can approximate the critical point 0(7) based on the asymptotic result 
in Theorem [TJ We can expect that this approach yields good approximation since in testing 
hypothesis (fTBT) . we often encounter the situation where the eigenvalues \ m+ i, . . . , X p are much 
smaller than the other eigenvalues. The hypothesis fll6p with small a 2 corresponds to the case 
im+i = • • • = £ p = 1 in Theorem [TJ Consequently we can approximate the distribution of V in 
(fT7|) by the distribution of 

TT P d- 

v= l p u=m+1 y . m (P-mr m , 

where dj (i = m + 1, . . . ,p) are the eigenvalues of Wishart matrix W p - m (n — m, J p _ m ). Even under 
the distribution Wp_ m (n — m, I p - m ), if is not easy to derive analytical expressions for percentage 
points for V. For p — m = 2, the distribution function is explicitly given (see 10.7.3. of Anderson 
(2003)) by 

P(y <v) = v^-™- 1 )/ 2 

which gives the critical point 0(7) as 

C ( 7 ) =7 2/(n-m-l) > ( lg) 

Generally a numerical calculation is needed for the exact evaluation of critical points. For this 
problem, refer to Consul (1967) and Pillai and Nagarsenker (1971). 

We made a simulation for the comparison between the above two methods. Let p = 3, m — 1 
and consider testing the hypothesis A2 = A3. We examined the accuracy of the two critical points 
(|18p and f[T9"j) with 7 = 0.05 and 7 = 0.01 by simulating the probability of V being smaller than 
these critical points when A2 = A3, that is, the probability of the error of the first kind. We put 
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Table 2: Simulated Type 1 Error 






1 


.9 


.8 




5%1 


5%2 


1%1 


1%2 


5%1 


5%2 


1%1 


1%2 


5%1 


5%2 


1%1 


1%2 


n = 5 


.140 


.041 


.051 


.008 


.142 


.041 


.053 


.008 


.141 


.042 


.054 


.008 


n = 10 


.063 


.033 


.015 


.005 


.064 


.033 


.016 


.006 


.067 


.036 


.017 


.007 


n = 20 


.038 


.026 


.007 


.004 


.039 


.028 


.008 


.005 


.041 


.030 


.008 


.005 


n = 50 


.025 


.021 


.004 


.003 


.027 


.024 


.004 


.003 


.032 


.029 


.006 


.005 


n = 100 


.022 


.020 


.003 


.003 


.025 


.023 


.003 


.003 


.034 


.031 


.005 


.005 


n = 500 


.016 


.016 


.002 


.002 


.029 


.029 


.005 


.005 


.043 


.043 


.008 


.008 


n = 1000 


.017 


.017 


.002 


.002 


.035 


.035 


.006 


.006 


.048 


.047 


.009 


.009 






.6 


.5 


.3 




5%1 


5%2 


1%1 


1%2 


5%1 


5%2 


1%1 


1%2 


5%1 


5%2 


1%1 


1%2 


n = 5 


.145 


.044 


.055 


.008 


.148 


.044 


.056 


.009 


.158 


.047 


.059 


.009 


n = 10 


.072 


.039 


.019 


.008 


.076 


.041 


.021 


.008 


.086 


.046 


.023 


.009 


n = 20 


.051 


.036 


.011 


.006 


.057 


.042 


.012 


.005 


.067 


.049 


.016 


.009 


n = 50 


.046 


.040 


.009 


.008 


.053 


.047 


.011 


.009 


.056 


.050 


.012 


.010 


n = 100 


.050 


.047 


.010 


.009 


.050 


.047 


.010 


.009 


.052 


.049 


.010 


.009 


n = 500 


.050 


.050 


.010 


.009 


.050 


.050 


.010 


.010 


.052 


.051 


.010 


.010 


n = 1000 


.051 


.050 


.010 


.010 


.051 


.050 


.010 


.010 


.051 


.051 


.010 


.010 







.1 


.01 


.001 




5%1 


5%2 


1%1 


1%2 


5%1 


5%2 


1%1 


1%2 


5%1 


5%2 


1%1 


1%2 


n = 5 


.161 


.049 


.061 


.010 


.164 


.050 


.063 


.011 


.167 


.051 


.064 


.010 


n = 10 


.091 


.051 


.026 


.011 


.092 


.051 


.025 


.009 


.090 


.050 


.024 


.010 


n = 20 


.066 


.049 


.015 


.010 


.067 


.050 


.016 


.010 


.067 


.050 


.016 


.010 


n = 50 


.057 


.051 


.012 


.010 


.057 


.051 


.012 


.010 


.055 


.048 


.012 


.010 


n = 100 


.053 


.050 


.011 


.010 


.053 


.049 


.011 


.010 


.054 


.051 


.011 


.010 


n = 500 


.051 


.051 


.010 


.010 


.050 


.049 


.010 


.010 


.051 


.050 


.011 


.010 


n = 1000 


.050 


.050 


.009 


.009 


.051 


.051 


.010 


.010 


.051 


.050 


.010 


.010 



Ai = 1, A2 = 0i A3 = (3 and varied both (3 and n. Table 2 shows the result, where the labels 
5(1)%1 and 5(1)%2 indicate that the numbers below correspond to the critical points ffTB"]) and ffl~9]) 
respectively with 7 = 0.05(0.01). Numbers in bold mean that they are within ±0.001 deviation 
from the desired value. 

We can summarize the result as follows; 

1. If > 0.8, both ffTB"]) and (fl~9l) need a large sample size. Especially when is as large as 
1.0 or 0.9, more than 1000 samples are required for a good approximation. There is no 
meaningful difference between both critical points. 

2. If equals 0.6 or 0.5, 50 (sometimes 20) samples are large enough to give a good approx- 
imation for both f fTB]) and ( jl~9l . There is no significant difference between both critical 
points. 
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3. If (3 < 0.5, f|T9|) shows significantly better performance than ([TBI . Even with such a small 
sample as 5, (fT9l gives very accurate approximations. The critical point (fT9l is robust in 
the sense that it already gives an excellent approximation when the smallest eigenvalues are 
0.3 times as large as the largest eigenvalue. 
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